Current correlations and quantum localization in 2D disordered systems with broken 

time-reversal invariance 
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We study long-range correlations of equilibrium current densities in a two-dimensional mesoscopic 
system with the time reversal invariance broken by a random or homogeneous magnetic field. Our 
result is universal, i.e. it does not depend on the type (random potential or random magnetic field) 
or correlation length of disorder. This contradicts recent tj-model calculations of Taras-Semchuk 
and Efetov (TS&E) for the current correlation function, as well as for the renormalization of the 
conductivity. We show explicitly that the new term in the cr-model derived by TS&E and claimed to 
lead to derealization does not exist. The error in the derivation of TS&E is traced to an incorrect 
ultraviolet regularization procedure violating current conservation and gauge invariance. 



The quantum coherence is well known to play a cen- 
tral role in the physics of mesoscopic systems. It in- 
duces, in particular, long-range spatial correlations of lo- 
cal densities of states and eigenfunction amplitudes. The 
long-range character of these correlations is due to the 
existence of massless modes, diffusons and Cooperons. 
This leads also to strong mesoscopic fluctuations of global 
quantities, such as the conductance or the inverse partic- 
ipation ratio. 

In this paper, we present one more example of such 
long-range correlations induced by the phase coherence 
of a system. Specifically, we will calculate the correlation 
function of local equilibrium current densities of electrons 
subject to a (possibly smooth) random potential (RP) or 
random magnetic field (RMF). Non-zero local currents 
exist in a system if the time reversal invariance is broken 
by a magnetic field, either uniform or spatially fluctuat- 
ing. We will demonstrate that the result is independent 
of the type of disorder and has a rather universal char- 
acter. 

This paper is to a large extent based on results ob- 
tained by two of us (A.D.M. and P.W.) in 1994 in col- 
laboration with the late Arkadii Aronov. At that time 
our research was motivated by the paper where it 
was claimed that the above current fluctuations lead to 
delocalization of electrons in random magnetic field, in 
contradiction to our finding in |^ . Our results remained 
however unpublished. In fact, the unfortunate error in 
the calculation of Ref. was quite obvious (these au- 
thors confused the correlation function of equilibrium 
currents with the current response function), and we did 
not see a necessity to devote an additional publication to 
it. Also, evaluation of the correlation function of local 
currents is quite straightforward and follows essentially 
the same lines as the calculation of a typical value of the 
global persistent current 

However, in the meantime, the result of Ref. [Q has 
been quoted in a number of publications as one of the- 
oretical proposals concerning the problem of localization 
of two-dimensional (2D) fermions in a random magnetic 
field. Furthermore, the idea of pi was recently revived. 



within a different line of reasoning, by Taras-Semchuk 
and Efetov (TS&E) [Q. Despite the criticism of their 
work in the comment Jq] , TS&E have made public a 
rather extended paper [^presenting in detail their calcu- 
lations and conclusions. If true, their results, predicting a 
delocalizing one-loop contribution to the /3-function gov- 
erning the scaling of conductance for systems with bro- 
ken time reversal invariance (unitary symmetry class), 
would change completely our understanding of localiza- 
tion properties of these systems. In particular, they 
would force one to reconsider the phase diagram of the 
quantum Hall effect. Most noteworthy, the findings of 
TS&E challenge the notion of universality for this class 
of systems: the new contribution to the /3-function de- 
rived by them has a non-universal coefficient /3tse (see 
below) depending on the type (RP or RMF) and the cor- 
relation length of the disorder. We thus believe that it is 
timely to present our results on current density correla- 
tions and our criticism of Refs. [||J^ without the length 
restrictions implied by the Comment format. 

We consider the correlation function of local current 
densities ja{j^, E) at equilibrium in 2D (we set H — 1 and 
omit the spin degree of freedom). 
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where Gf, are retarded and advanced Green's func- 
tions, and (. . .) denotes disorder (RP or RMF) averaging. 
In view of the current conservation, Vj = 0, the function 
/a/3 is transverse in momentum space. 
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Since we are interested in long-range current correlations, 
we will study the low-momentum, low-frequency behav- 
ior of Ia/3. A finite value of F(q,Lo) in the low-g limit 
would imply dipole-type long-range correlations in coor- 
dinate space. 
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We will make use of the representation of the correla- 
tion function I^p in terms of a functional derivative of a 
generating functional ^{A,A'} with respect to vector- 
potential source fields A(r), A'(r), 



/a/3(r-r') 
where 



, (5) 



A=A'=0 



^{ A, A'} = ( [Tr In G^{A} - Tr In Gi{A}] 

X [TrlnGi+JA'}-TrlnG^+^{A'}]). (6) 

It is seen from Eq. (||) that the transverse character of 
/q/3 is intimately related to gauge invariance of the gen- 
erating functional J?^{A, A'}. 

We start from a diagrammatic analysis in the frame- 
work of a conventional impurity diagram technique. The 
diagrams for J-' are obtained by connecting two closed 
electron loops by impurity lines; the diagrams for I^p 
are then generated by inserting two current vertices (one 
in each loop) in all possible ways. 

Let us first consider the case of a random scalar poten- 
tial with a correlation function {U{r)U{r')) — M^(r — r'); 
its correlation length d is assumed to be finite, though 
it can be arbitrary as compared to the wave length Ap. 
In particular, if d ^ Af (smooth RP), the scattering is 
of a small-angle nature, implying that the transport re- 
laxation time is much longer than the single-particle one, 
Tti ^ It is clear that in any finite order n of the 
perturbative expansion [n is a number of impurity lines 
connecting the fcrmion loops), the current correlation 

fmiction is finite-ranged, i.e., it has a finite corre- 
lation length, beyond which the correlations decay expo- 
nentially in view of the finite-range character of W{r—r'). 

This means that in momentum space ij^^ (q) should have 
no singularity as g — > 0, implying that F{q) oc q^. In- 
deed, the sum of all the diagrams of the n-th order can 
be presented in the form 

iSlicO = ^ y (rfqi) ■ ■ • {dqn-l)W{qi) . . . W{q^) 

X T«(qi,...,q„)r0(-qi,...,-q„), (7) 

where (dq) = (Pq/ {2'it)'^ , and TQ(qi, . . . , q„) denotes a 
vertex part to which n impurity lines with momenta 
qi, . . . , q„ satisfying qi = q are attached in all pos- 
sible ways. Because of its vector character and the sym- 
metry with respect to permutation of q^, the block Tq 
should be proportional to q; the current conservation im- 
plies its transverse character oc Qq,, where Qq, = eafsQp 
[cafs is the antisymmetric tensor). When substituted in 
Eq. (^, this yields again 



We now turn to the RMF case. The situation now is 
somewhat less trivial, since although the RMF correla- 
tion function W_B(r — r') — {B{r)B{r')) is assumed to 
be of a finite range, the correlation function of the cor- 
responding vector potential A{r) is long-ranged, which is 
reflected by a singularity at g = in momentum space, 
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Let us, however, analyze the RMF analog of Eq. (|^), 

^5(q) = ^/(rfqi)---(rfq«-i) 

X WaiPi (qi) . . . Wa„p„ (qn) 
^ '^aa.i ...a„(Ql!) ■ ■ ■ ; Qn) 

X T^0i.../3„(-qi,...,-q„) . (10) 

As in the RP case, the block T represents a sum of all 
possible diagrams with one current vertex and n disor- 
der lines attached. This includes, in addition to the di- 
agrams appearing in the RP case, diagrams having ver- 
tices originating from the A^ term in the Hamiltonian, 
H = (-iV - eA(r)/c)2/2TO (i.e. with two RMF lines 
joining the electron line at the same point, or with a 
RMF line starting at the external current vertex). Such 
diagrams, while being of minor importance for the diffu- 
sion contribution considered below (they are suppressed 
by the factor (fcp/tr)"^, where fcp is the Fermi momentum 
and /tr = VFTtrYj are crucially important in the pertur- 
bation theory 0. It is not difficult to see that all the 
diagrams for the vertex part T in Eq. ( p^ ) can be gener- 
ated by variation of a current with respect to the vector 
potential 



^aai ...an (qi 5 ■ ■ ■ i qn) 
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A=0 



In view of the gauge invariance, the vertex part in the 
limit q ^ must have the form 

{<ll,---,<ln) 0Cqa{qi)ai---{q7i)an ; (12) 

SO that all the singularities in the vector potential correla- 
tors WaiPiiqi) are canceled by the vertex parts, yielding 
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FRp{q) oc q^ 
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Therefore, in any order n the perturbative contribution 
to the current correlator is of finite range in the RMF 
case as well, despite the singular nature of the vector 
potential correlator (H). 

As usual in mesoscopic physics, long-range correlations 
are determined by a diffuson contribution. We assume 
Cooperons to be suppressed by the random or homoge- 
neous magnetic field; in the latter case the magnetic field 
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is assumed to be classically weak {wcTti <C 1). The cor- 
responding diagrams are shown in Fig. |^. The sum of 
these two diagrams satisfies the current conservation re- 
quirement, since they can be obtained from Eq. (||) with 
the generating functional J- shown in Fig. 0. 





FIG. 1. Leading diffuson diagrams for the current corre- 
lator with one (a) and two (b) diffusons. 



The diagram Fig. ^ gives 



^'/3 = / (dk)(dk')(dQ) 



X y„(k,Q,q)[l/;5(k',-Q,~q)]' 

X £'kk'(Q,t^)-Dkk'(Q + q-t^)> 

where Va(p) are the triangle vertex parts, 



(16) 



K,(k, Q, q) = t;„G«(k + q)G^(k - Q)G«(k). (17) 



Substituting (|lj) in (|16D, we find the following two rele- 
vant contributions. If one takes into account the singular 
parts of both diffusons, then the Green functions can be 
expanded in small momenta q and Q, which yields for 




FIG. 2. Generating functional T for the diagrams of Fig. |l| 

In order to take into account anisotropic scattering we 
use the following expression for the diffusons B 



1 7k(q)7k'(q) 



Ac 



(14) 



where = m/27r is the density of states, D ~ VpTti/2 is 
the diffusion constant, 7k (q) = 1 ^ ^(''"tr — T)kq/m, and 



Dice 



-x™(k)x;,(k'), 



with Xm(k), Am being the eigenfunctions and eigen- 
values of the scattering operator with the kernel 
VF(k, k')|G^(fc')|^. When calculating the contribution of 
the diagram Fig. |l|a, only the singular part of the diffu- 
son (the first term in (^4|) ) should be taken into account, 
yielding 



- / (dk)(dQ)w„«^[G^^(fc)G^(fc)]^Z?kk(Q,c^) 



(dQ) 
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where — {D/u})^/'^. Let us note that we consider an 
infinite system; for a finite sample of size L < one 
should replace by the system size L, which plays then 
the role of the infrared (IR) cut-off. 



27r2 



idQ){2Q^ + q^){2Qp + qp) 
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Another contribution 
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(18) 



(19) 



arises when the regular part of one of the two diffusons 
is taken into account. 

Combining Eqs. ([s]), (|l|) and we finally get 



Lap 



■ ln(gic. 



'a/3 



9a 9/3 



27r3 V r 

which has the required transverse form 



(20) 



with 



L-J « q « l^,^ . (21) 



Note, that neither the single-particle nor the transport 
scattering time enter the final result, which is universal 
and does not depend on the character of the disorder. 

Our result (|o|) differs from that of TS&E (Ref. §) in 
two crucial aspects. First of all, in the formula by TS&E 
the factor F{q,ijj) ^ h\{qL^^ is replaced by a large con- 
stant /3tse > 1- Specifically, for the RMF case TS&E 
find 



R Ttr 



(22) 



while for a smooth RP their consideration yields Q 

/3tse^ (^)'^'. (23) 



Secondly, while our contribution comes from the IR re- 
gion Q <C l'^^ ^ their term is claimed to be determined by 
the ultraviolet (UV) region Q ^ V^^ . We will now show 
that the result of TS&E is wrong. 



3 



In fact, already the above considerations imply that the 
findings of TS&E are incorrect. Indeed, we have demon- 
strated that F{q,uj) oc as (7 ^ In other words, 
Eq. (|^) implies F{q — 0,uj) — 0, since a gauge-invariant 
generating functional cannot depend on a change of 
the vector potential by a constant. Therefore the TS&E 
result F{q = 0,0;) = /3tse ^ 1 cannot be correct. 

To make closer contact with the work of TS&E and, in 
particular, to show explicitly their mistake, we will now 
rederive Eq. ( po|) i n the framework of the ballistic a- 
model approach |ll| , p^ . The cr- model can be derived by 
averaging over a smooth RP or a RMF following |^,|| , as 
explained in detail in . Performing a standard calcula- 
tion and treating the cr-model correlator perturbatively, 
one gets the leading contribution in the form 



Iai3{r - r',uj) = 



dndn' 
2^2 "^F / j^:;^ 



i2ny 

xn„5(r, n; r', n'; uj)n[3g{r' , n'; r, n; w), (24) 

where ^(r, n; r', n'; lo) is a propagator describing the mo- 
tion from the point (r', n') of the phase space to the point 
(r, n), with the unit vector n staying for the velocity di- 
rection of a particle on the Fermi surface. The propagator 
g obeys the Liouville-Boltzmann equation 



(£ - iuj)g{r, n; r',n';uj) 
27r 

5(r - r')(5(n - n') 



'"(n, n")[g(r, n"; r', n'; uj) - g{r, n; r', n'; u)] 

(25) 



where C = vpiiVr and 'w{n,n') is the scattering cross- 
section {(j) is the angle between n and n'). 



w{n, n') 



27rz^VK(2pF sinf), 
27riycot2 | 

2 V 2mc ) 



iyB(2pF sin- 



RP, 
RMF. 

(26) 



Note that (gj) is identical to the Altshuler-Shklovskii- 
type contribution to the local density of states correlation 
function |p2[, up to the velocity factors v-pnonV-pUfi. 

Equation (|2J) has the same form as the formula (37) 
of TS&E, Rcf. 1^. In fact, there is a slight difference 
between the two formulas: while our function 5 is a 
full propagator of the Liouville-Boltzmann equation, the 
function r(r, n; r', n') of TS&E is obtained by projecting 
this equation onto the space of non-zero harmonics in the 
velocity space. As a consequence, our propagator has a 
diffusive behavior at large distances (see below), while F 
of TS&E does not. This is because we treat the cr-model 
field fully perturbatively, while TS&E do this with re- 
spect to the non-zero harmonics only. Therefore some of 
the contributions to our Eq. (|2j) are shifted by them into 
other terms not included in Eq. (37) of Ref. While 
being crucial for a proper account of the diffusive contri- 
bution, this difference is irrelevant for the analysis of the 
contribution of the ballistic region, with characteristic 



momenta of the propagators satisfying Q ^ l^^ , which 
is claimed by TS&E to be responsible for the results of 
Ref. |§. 

In the hydrodynamic limit (QZtr ^ 1, i^^tr ^1) the 
propagator (7(Q; n, n'; w) has the conventional diffusive 
behavior (in view of the vector nature of the factors 
UajUf} in ( p^ ) we have to take into account the leading 
correction in Qhr), 

ff. , ^ (l-^^trQn)(l-Z^trQnO 

gQ;n,n;t^= — ■. . 27 

DQ'^ — luj 

Substituting ( p7|) in (p4|), we arrive again at Eq. (|20|). 

Let us now demonstrate how the statement F(g — 
0, w) = is obtained within this approach. This can 
be most clearly done if one introduces the path integral 
representation for the propagator 

g(r,n;r',n';c^)= / dTe^"^g(r, n; r', n'; T); 
Jo 

.z(T) = {r,n} 

g(r,n;r',n';T)= / Vv{t)e^'^<''>\ (28) 

Jz(0) = {r',n'} 

where r(i) = UFn(t) and z(i) = {r(i),n(t)}. Here S is 
the action corresponding to the dynamics determined by 
Eq. (Eq). Using (Eq), it is easy to show that 



y Av'g{v, n; r', n'; tj)n^g(r', n'; f , n; w) 



00 <.z'(T) = {r.n} fT 

dTe'^^ / Vr'{t) / dtn'^{t)e'^^''''-'^'i 

Jz'(0) = {f,n} Jo 



.g(r,n;f,n;w; 



VF 



(29) 



Substituting this (with r = f ) in Eq. ( p4[ ) integrated over 
r', we immediately get lapil = 0,w) = 0, in agreement 
with our earlier analysis. This can be traced back to the 
fact that an integral of velocity vector over a closed loop 
is equal to zero, see Fig. ||a. Note that one should exclude 
from Eq. ( p9| ) a contribution of a straight trajectory of 
zero length in the limit f ^ r (see a detailed discussion 
below). In fact, this proof of F{q = 0, cj) = can be cast 
in the form fully analogous to the above diagrammatic 
derivation based on Eq. (|). Indeed, Eq. (|^) can be 
rewritten in the form (|5|) with the generating functional 
given by an integral over closed paths in the phase space, 



^{A,A'}=/ 



z(0)=z(T) 



Vv{t) 



:exp [iS[r{t)] ~ i-^ j dtr{t) [A(r(t)) + A'(r(t))] 



(30) 



The trajectories rit) in Eqs. (|29|), ( pO[ ) are classical coun- 
terparts of the diffuson ladders in Fig. |l]. Since the in- 
tegral in (^) goes over closed loops, the functional T 



4 



is gauge-invariant (the second term in the action yields 
the flux through the loop), which immediately leads to 
F{q — 0, Lu) — as explained above. 





(b) 

FIG. 3. (a) The paths contributing to the current corre- 
lation function I^p , Eq. (p^; (b) paths giving rise to the 
spurious contribution to 1^^ found by TS&E. 

Why do TS&E fail to observe that F{q = 0, lu) is iden- 
tically zero and find instead the results (p2[), (23) for 



this quantity? The reason is as follows. TS&E evaluate 
Eq. ( p^ by writing it in momentum space and imposing 
an UV cut-off Qq ^ 1^^.^ which is set Qq = (upr)"^ in 
the end of the calculation. It is easy to see that within 
this procedure Eq. (p3) is replaced by 



-'a/3 \^ 



27r2 



df A(r - f , Qo)Ja(3ir, r', f ; uj), 

(31) 



where 



Ja/3(r,r',f;cj) = J 



dndn' 



(2^)2 

XTiagir, n; r', n'; uj)n'pg{r' , n'; f , n; uj) 
and A(p, Qq) is a broadened (5-function, 



A(p,Qo)= / (dQ) 

J\Q\<Qa 



QpJiiQop) 

2'Kp 



(32) 



(33) 



with Ji (x) the Bessel function. The crucial difference be- 
tween Eq. (|3l| ) and the exact formula (|2j) is that the for- 
mer is determined by non-c/os erf paths withjf — r| ^ Qn^ ■ 
In particular, there is a contribution to (^ij) from a very 
short, almost straight path of a length ~ |f — r| , with the 
point r' lying between the points r and f (see Fig. ||b). 
It is precisely this spurious contribution that determines 
the result of TS&E. To demonstrate this, we simply cal- 
culate this contribution. According to (|3l|), evaluation of 
the q ^ limit of I^^^iQ^ ^) involves integration of Jap 



over r'. Using (|29|), we easily find that the contribution 
of the almost-direct trajectories to J dr' J^p is equal to 



j dr'Jai3{r,r',r;u;)=J ^^(r, n; f , n; w) 



PaPf3 

vfP 



, (34) 



where p ~ r r. If one completely neglects the disor- 
der, g becomes the free propagator (since we discuss a 
ballistic contribution, the frequency w is fully irrelevant 
assuming u; < t^^ , and one can simply set w = 0), 



ffo(r, n; r',n') 



1 



27rwF|r — I"' I 



tS n 



r - r' 



5{n-n'), 
(35) 



substitution of which into ( p4|) leads to a divergency of 
the type 6(0). To see how this divergency is smeared by 
the disorder scattering, we rewrite (p4f) in the form 



/ dr'XMr,r',f;c.) = Pf§r{0, p/vp), 
J v^P 



(36) 



where 'P{S(j),t) is the distribution of the variation dcf) of 
the velocity angle in time t. In the case of a RP the 
velocity angle diffuses with the diffusion constant r^^^, 



■PRp{6^,t) 



exp 



" 4t 



while in a RMF one has a Cauchy distribution ||l^] 
1 2i/rtr 



VRMF[54i,t) 

Substituting 



TT (2t/Tt,)^ + (-50)2 ■ 

81) in (|3^) we arrive at 



(37) 



(38) 



J dr'Jai3ir,r' ,r;Lj) 



1 PcPfS ^ / (^t,■/4^p)^/^ RP 



vl p2 I k,/2TTp, RMF. 
Finally, substituting (pS) in (pH), we get 



(39) 



7(Q0^tr)^/', RP 

Qo^tr/4, RMF, 



(40) 



with a numerical coefficient 7 — ■\/7r/2[r(|)/4r(|)], re- 
producing the resuhs (H), ^) of TS&E. 

We have thus demonstrated that the findings of TS&E 
result from an incorrect regularization procedure. Note 
that their term ( ^ ) is not transverse with respect to q 
and thus violates the current conservation. One could an- 
ticipate this without any calculations: a non-zero q — > 
limit of in combination with current conservation 
would imply, according to (^, existence of a singular 
term ^ qalp/l^ of a long-range nature (^ in coordinate 
space. However, TS&E deal with short scales ^ hr only 
and thus have no chance to obtain such a term. This 
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violation of current conservation remains, however, hid- 
den in 1^, where only the trace laa^iQ = 0,a; = 0) is 
calculated, while the transverse structure of is then 
"restored" by hand. 

The violation of current conservation is a direct conse- 
quence of the violation of gauge invariance of the generat- 
ing functional by the regularization of TS&E. Indeed, 
this gauge invariance is guaranteed by the fact that only 
closed paths contribute to (|3^) while TS&E allow for non- 
closed short paths which produce their spurious term. It 
is clear from our discussion what would be a proper UV 
regularization procedure respecting the gauge invariance 
constraint. One should discard the contribution of those 
closed loops to the effective action which have a length 
shorter than a cut-off length vpTg^ which amounts to 
setting the lower limit of the T-integral in (^o|) equal to 
Tq. Note that for a smooth disorder {d 3> Xp) the qua- 
siclassical result ( p^ ) for the current correlator will be 
essentially independent of the cut-off To ^ rtr, since a 
particle cannot make a closed loop in a time much shorter 
than Ttr (if one neglects exponentially rare events). 

To shed more light on the cause for the appearance of 
the contribution /J^^, let us note that the length of the 
relevant trajectories p ^ 1/Qo goes to zero when the cut- 
off is removed, Qq oo. Such "zero-length trajectories" 
are well known to appear in the quasiclassical represen- 
tation of the Green's function G£;(r, f) at f ^ r. It is 
also well known that the contribution of such a direct 
trajectory of length zero should be split off from the re- 
maining sum over periodic orbits in the trace formula for 
G_E(r, r); the former yields the average (Weyl) contribu- 
tion to the density of states, while the latter describes 
fluctuations around it |l^] . In the present context of cur- 
rent correlations, the contribution of the zero-length tra- 
jectory is zero since (j{r,E)) = (it corresponds to the 
disconnected part of lajs described by a diagram with no 
impurity lines at all). The rest is represented by a sum 
over closed classical paths (Eq. ( p^ ) and Fig. ||a) with the 
direct zero-length trajectory excluded. This is in full anal- 
ogy with the semiclassical consideration of the spectral 
correlator {v{E)v{E+lo)) jl^, where the zero length tra- 
jectory yields the disconnected part {v)"^ , while the sum 
over periodic orbits reproduces the Altshuler-Shklovskii 
diffuson contribution. Note that the procedure used by 
TS&E would produce an additional spurious contribu- 
tion of very short, almost straight paths in this case as 
well. 

We now turn to the closely related analysis (already 
presented in a brief form in the Comment ^) of the new 
term in the cr~model obtained by TS&E in Rcfs. [^,1). 
According to these authors, this term leads to delocal- 
ization at the one loop level, so that the corresponding 
diagrams should contain only one true (singular at low 
momenta and frequency) diffuson. One can show, how- 
ever, that, very generally, such a contribution to con- 
ductivity does not diverge. For white-noise RP this was 
shown in Ref. M]. We will now prove that the state- 



ment holds also in the case of RMF or smooth RP with 
time reversal invariance broken. Indeed, the sum of all 
one-diffuson diagrams with n impurity lines crossing the 
diffuson and having their starting (end) point anywhere 
on the left (right) block has the form 



27r 



(dq) 



Mi"'(q,c^) 



2'KVT'^{Dq'^ — iuj) ' 



(41) 



where 



(dqi)...(dq„_i)I^((Zi)...T4^(q„) 



X Sl{cii,...,q,n)Sf{-q,i,...,~Cir,), RP 
M(S)(q,c.) = l y'(dqi)...(dq„_i) 

X >VQi/3i(qi) • ■ ■ >Va„/3„(qn) 
X 'S'xai...a„(qi,---:qn) 

X ^53i.../3„(-qi'---'-q«)' ™- (42) 

Here 5''^ [S^) denotes the left (right) vertex part which 
has one external current vertex and to which n impurity 
lines (momenta q^) as well as one diffuson (momentum 
q = ^ are attached. Using the same arguments as 
presented above for the vertex part T of the current corre- 
lator, one finds |] S^-^ cx q for RP and S ' oc qqi . . . qn 
for RMF, so that the correction ( ^l|) to conductivity is 
non-divergent as w — > 0, independently of the type of dis- 
order and presence or absence of the time-reversal sym- 
metry, 



A4:)(^)a (dq) 



{Dq^ - iuj) 



cx const - 



(43) 



In fact, the absence of a divergent correction is inti- 
mately related to the above analysis of the q ^ form 
of the current correlator. Indeed, summing up diagrams 
with a diffuson vertex inserted in all possible ways and 



usmg 



finds 



e'Mi^^i>{q = 0,c.) = {27TTflil>{q = 0,c.), (44) 

so that the prefactor in front of the logarithmic diver- 
gency ^ {dq) / v{Dq^ ~iuj) in (^) is equal to llcx (q — 0,ct;) 

(as also found by TS&E). Since, however, ll,^ [q — 
0, uj) = 0, this just amounts to saying that the diver- 
gent one-diffuson contribution does not exist, neither for 
the RP nor for the RMF problem. As to the error in the 
cr-model calculation of conductivity by TS&E in |^,^ , it 
is identically the same as in their calculation of the cur- 
rent correlator (discussed above in great detail), in view 
of the relation (H) between the two quantities. 

As is well known, divergent contributions to the con- 
ductivity appear if one allows for more than one dif- 
fuson. In particular, the counterpart of the one- and 
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1 



(45) 



two-difFuson diagrams for the current correlator (Fig. |^) 
are the conventional two-loop diagrams for conductivity 
with two and three difFusons, yielding the usual weak- 
localization correction in the unitary ensemble, 

A5 _ 1 
9 

where g = a^x/ {e-'^ /h). 

Finally, let us return to the one-diffuson corrections 
to conductivity (^). Though non-divergent, they are 
still of considerable interest, because of their non- analytic 
(oc \uj\) character at w ^ 0. Such contributions corre- 
spond to \/t^ long-time tails in the velocity correlation 
function and are induced by return processes neglected 
in the Boltzmann transport theory. They were discov- 
ered long ago in the context of the classical Lorentz gas 
10 and were recently studied for a 2D electron gas 
subject to a smooth RP or RMF. It was found that, in 
contrast to the universal weak localization corrections, 
the prefactor and the sign of the |a;| term is non- universal 
(i.e. depends on the type of disorder). 

In conclusion, we have studied long-range correlations 
of local current densities j(r, E) in a 2D mesoscopic sys- 
tem with time reversal invariance broken by a homoge- 
neous or random magnetic field. We have considered 
two types of disorder: a (possibly smooth) random po- 
tential and a random magnetic field. We have demon- 
strated that within any given order of the perturbation 
theory the range of correlations remains finite even in the 
RMF case, despite the IR-singular nature of the vector 
potential correlator. The long-range correlations are de- 
termined by the diagrams with difFusons, Fig. |^, yielding 
the result (pO|). The current correlation function is found 
to be universal, i.e. it does not depend on the type or 
correlation length of disorder. 

Our results demonstrate that recent findings of TS&E 
1^,^ on the current correlations and on the conductivity 
renormalization leading to delocalization in 2D electron 
systems with a smooth RP or a RMF (as well as a sim- 
ilar statement of an earlier paper Q ) are fundamentally 
wrong. To make this point as clear as possible, we have 
demonstrated how our results are reproduced within the 
ballistic cr-model approach. We have further shown that 
the spurious contribution found by TS&E results from 
an incorrect regularization procedure violating current 
conservation and gauge invariance. This term does not 
arise if one uses diagrammatic perturbation theory or a 
properly regularized ballistic cr-model. Therefore, the 
new term in the diffusive cr-model which was obtained 
in and which was claimed there to lead to delocal- 

ization, does not exist. In particular, the RMF problem 
belongs to the conventional unitary symmetry class with 
the leading quantum correction being of two-loop order 
and of localizing nature as found earlier . 
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